STEADY WATER WAVES WITH MULTIPLE CRITICAL LAYERS 



MATS EHRNSTROM*, JOACHIM ESCHER+ , AND ERIK WAHLENt 

Abstract. We construct small-amplitude periodic water waves with multiple critical layers. In 
addition to waves with arbitrarily many critical layers and a single crest in each period, multimodal 
waves with several crests and troughs in each period are found. The setting is that of steady two- 
dimensional finite-depth gravity water waves with vorticity. 

1. Introduction. The periodic steady water-wave problem describes wave-trains 
of two-dimensional surface water waves propagating with constant shape and velocity. 
The mathematical existence of such steady water waves is well-established, and one 
distinguishes between infinite and finite depth, gravitational and capillary waves, and 
flows with and without vorticity. Historically, the main focus has been on irrotational 
waves, i.e. on flows without vorticity. An important class of such waves is given by the 
Stokes waves: periodic gravity- waves with a single crest in each period around which 
the wave is symmetric. Of particular mathematical interest is the so-called wave of 
greatest height with a sharp crest of interior angle 120 degrees. The paper [21] de- 
scribes much of that theory; see also [14] for a general overview of the water-wave 
problem. 

A few years ago the authors of [7] proved the existence of water waves with a 
general vorticity distribution (for a different approach, cf. [20]). This opened the 
way, not only for a new branch of research, but also for the existence of new types of 
waves. In many aspects, the features of rotational waves resemble those of irrotational 
waves: results on symmetry [4, 5, 15], regularity [6], and waves of greatest height [23] 
have now been extended from irrotational waves to waves with a general vorticity 
distribution. In most of those cases the differences have manifested themselves as 
technical difficulties, not different properties per se. 

An exclusion from this rule is the appearance of internal stagnation points — 
fluid particles with zero velocity relative to the steady flow — and critical layers. By 
a critical layer we mean a connected part of the fluid domain consisting only of 
closed streamlines. The existence of such waves within linear theory was established 
in [13], where the streamlines and particle paths of waves with constant vorticity were 
studied. For irrotational flows such behavior is not encountered, cf. [3, 8]. Also earlier 
studies had indicated that internal stagnation should be possible (see, e.g., [16]), but 
it was first with [26] that waves with critical layers were given a firm mathematical 
foundation (see also [9]). A different approach in dealing with the free boundary 
allowed for including internal stagnation, a phenomenon which for technical reasons 
had not been possible to see with earlier methods for rotational waves (cf., e.g., [7, 24]). 

Another feature that separates waves with vorticity from those without is that 
the corresponding linear problem allows for multi-dimensional kernels. Whereas cap- 
illarity naturally introduces more degrees of freedom in the linear problem [18, 17, 
22, 24, 25, 27], periodic irrotational gravity- waves allow only for kernels with a single 
Fourier mode [12]. Using a result from inverse spectral theory, it was shown in [12] 
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that one can construct rotational background flows for which the kernel of the linear 
problem contains several Fourier modes, even in the absence of surface tension. It 
was, however, not clear to what extent those flows would be "natural" or "exotic", or 
whether they would allow for exact nonlinear waves. 

In this paper we extend the results from [26] and combine them with those from 
[12]. The novelty is two- fold: we construct i) waves with arbitrarily many critical 
layers, and ii) families of waves bifurcating from kernels with two different Fourier 
modes. The latter type of wave displays a more complicated surface profile than 
those of Stokes waves (cf. Figure 5.1, p. 19). This qualitatively different behavior 
arises not as a result of involved rotational flow, but the vorticity functions are linear 
maps. Thus, the paper at hand is a natural extension of [26], which dealt with flows of 
constant vorticity. We settle here for the mere existence of such waves. A qualitative 
investigation of the small-amplitude waves corresponding to the first Fourier mode 
will be given in a separate study [11]. 

For a given basic wave number k > 0, our problem has three non-dimensional 
parameters. Treating the first of these as a bifurcation parameter, while fixing the 
other two, we prove that a curve of solutions of the nonlinear problem bifurcates from 
any point where the kernel of the linearized problem is one-dimensional using the 
bifurcation theory of Crandall and Rabinowitz [10]. 

At points where the kernel is two-dimensional we similarly fix the third parameter 
and treat the other two as bifurcation parameters. Assuming that the linearized 
problem has two solutions with wave numbers fci , &2 G kZ + , such that k\ < k^ and 
kilk\ ^ Z, we obtain in this way three two-dimensional sheets of bifurcating solutions. 
The first and the second of these two sheets consist of "regular" solutions of minimal 
period 2n/ki and 2-7r/fc2, respectively, with one crest and one trough per period, while 
the third consists of multimodal solutions, that is, solutions with several crests and 
troughs per period (and with minimal period greater than 2ir/ki), except for the 
points where the sheets intersect. The points of intersection can be interpreted as 
secondary symmetry-breaking bifurcation points along branches of regular solutions 
with period 2n/ki or 27r/fc2- We also obtain some results in the more complicated 
case when kijk\ <E Z. Our approach is based on the Lyapunov-Schmidt reduction 
and the implicit function theorem. The same method was recently used to study the 
bifurcation of steady irrotational water waves under a heavy elastic membrane [1]. 

For a comparison with the irrotational case, we would like to draw the attention 
to the investigation [2] . There a global bifurcation theory for Stokes waves was devel- 
oped. The authors prove that the main bifurcation branch of, say, 27r-periodic Stokes 
waves contains countably many subharmonic bifurcation points, from which waves 
with minimal period 27rg, q a prime number, bifurcate. These bifurcations do not 
take place in the vicinity of laminar flows, but as one approaches the wave of greatest 
height along the main bifurcation branch. In the presence of surface tension, how- 
ever, the linear problem has a two-dimensional kernel for certain critical values of the 
surface tension parameter. The associated bifurcation problem has been studied in 
detail [18, 17, 22]. Using the Lyapunov-Schmidt reduction and blowing-up techniques 
the problem can be reduced to a finite-dimensional system of polynomial equations. 
For fixed values of the surface tension, branches of regular solutions bifurcate when 
the wave speed is used as bifurcation parameter. Multimodal solutions appear ei- 
ther as primary bifurcation branches or through secondary bifurcations from primary 
branches of regular solutions, depending on the value of the bifurcation parameter. 

The plan of this paper is as follows. Section 2 describes the mathematical problem, 
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the laminar background flows and the transformation of the free-boundary problem 
to a fixed domain. In Section 3 we fix the functional-analytical setting and introduce 
some auxiliary maps. Section 4 contains the main results, which in Section 5 are 
applied to yield three different families of solutions to the original problem. 

2. Preliminaries. The steady water-wave problem is to find a surface profile r\ 
and a velocity field (u,v), defined in the fluid domain = {(x,y) G R 2 : < y < 
d + Tj(x)}, where d is the depth of the undisturbed fluid, satisfying Euler's equations 

(u - c)u x + VUy = -P x 

mfl, (2.1a) 

(U - c)v x + vv y = -P y - g 

and 

u x + v y = mil, (2. lb) 

together with the kinematic boundary conditions 

v = at y = (2.1c) 

and 

v = (u — c)rj x at y = d+rj(x), (2-ld) 

as well as the dynamic boundary condition 

P = at y = d + rj(x). (2.1e) 

Here P is the pressure, c > is the wave speed and g is the gravitational constant of 
acceleration. 



Stream-function formulation. Denote by B the flat bed y = and the by S 

the free surface y = d + r)(x). The continuity equation (2.1b) implies the existence 
of a stream function if) with the properties that ip x — —v and ip v = u — c. The 
vorticity of the flow, uj :— v x ~ u y , describes the rotation of the water particles. Note 
that u> = —ipxx — ipyy The assumption that u < c in Q implies the existence of a 
function 7, such that ui = j(tp) (see [7]). Under this assumption, the equations (2.1) 
are equivalent to the following boundary value problem for ip: 

4>xx + Ipyy = ~7(V0 

i> = mo 
%j) = mi 

with Tp y < in £1. Here too and mi are real constants and the vorticity function 7 
is an arbitrary real- valued function defined on the interval [mi , too] . The difference 
mi — too is the relative mass flux and Q is a constant related to the energy. Dropping 
the assumption that u — c < one still finds that a solution of the problem (2.2) 
generates a solution of (2.1). The simplest example of a non-zero vorticity function 
is a constant function. This case was first investigated in [13] for linear waves, and 
then put on a rigorous foundation in [26]. Here we treat the second easiest example, 
namely that of an affine function 7. This choice of 7 is primarily motivated by its 
mathematical simplicity rather than by any physical arguments. However, an affine 
function can of course be seen as a Taylor approximation of a more general 7 and in 
this sense our choice seems natural. 
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in f2, 

on B, 

on S, 

on 5, 



(2.2) 



Linear vorticity function. Notice that (2.2) with an affine, non-constant vor- 
ticity function can be reduced to the same problem with a linear vorticity function 
by redefining ip. We therefore consider the following problem, where a is an arbitrary 
constant and we have normalized gravity to unit size: 

tpxx + ipyy = otip in £7, 

ip = m n on B, 

, a (2-3) 

ip = mi on 



1 IV7„/.|2 
2 



|Wr+r? = <2 on S. 



We also normalize the undisturbed depth to unity, so that the surface S is given by 
y = 1 + rj(x) and the fluid domain by £1 = {(x, y) € R 2 : < y < 1 + ?l(x)}. We 
pose the problem for n e C 2 +^(§, K) and ip € Cp+^ ovon (f2, K), where the subscripts 
per and even denotes 27r/K-periodicity and evenness in the horizontal variable, and 
we have identified 2-7r/K-periodic functions with functions defined on the unit circle S. 
Here, k is a positive real number. We also require that min?? > — 1. The parameter a 
influences the nature of (2.3) in a substantial way. For example, the trivial solutions 
are monotone in y when a > whereas they oscillate when a is sufficiently negative 
(cf. (2.4)). In order to simplify the discussion below we therefore make the following 
assumption. 

Assumption 2.1. We assume that a < 0. 
When a = the vorticity function is constant and one cannot transform (2.2) to an 
equivalent problem with a linear vorticity function. This case is, however, covered 
in [26]. When a > there are only simple bifurcation points (that is, with one- 
dimensional kernels). It is not difficult to prove that a curve of nontrivial, regular 
solutions bifurcates from each of these points, just like in Theorem 4.6. 

Laminar flows. The solutions constructed in this paper are perturbations of 
laminar flows in which the velocity field is horizontal but depth-dependent. The 
laminar flows are trivial solutions of problem (2.3) in the sense that rj = and ip is 
independent of x. Mathematically, they are given by the formula 

ipo{y,iJ,,X,a) := fj,cos(6 (y- 1) + A), (2.4) 

where /i,Ael are arbitrary constants and 

k ■= |a + fc 2 | 1/2 , keR. 

The numbers Q = Q(fj,,a, A), mo = mo(p, a, A) and mi = mi(p, X) are determined 
so that the function ip defined by (2.4) solves equation (2.3), i.e. 

/i 2 flo sin 2 (A) 

Q{Hi A, a) := (2.5) 

and 



mi(/i, A) :~ ficos(X), mo(/i, a, A) := /icos(A — #o) 

We make the following assumption on the parameters. 

Assumption 2.2. We assume that (ipo) y (l) ^ 0, that is //6> sin(A) ^ 0. 
It can be shown that without this assumption the linearized operator D w T appearing 
in the bifurcation problem is not Fredholm; see Remark 3.7. 
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The flattening transform. Problem (2.3) is a free-boundary problem in the 
sense that rj is a priori unknown. In order to find a secure functional-analytic setting 
in which to study the problem we use the change of variables 

(x,y) >-> (x 



1 + rj{x) 



which maps Q, onto the strip Cl = {(x, s) G R 2 : s <G (0, 1)}. In the coordinates (x, s) 
the problem (2.3) takes the form 

^-TT?J/iT^^-TT^J/(iT^ = ^ mn ' 

iJj = mn on s = 0, 

; (2.6) 

ip = mi on s = 1, 

+ 2TTT^ + r/ = on ' s = 1 - 

Let £((r),ip),a,Q) = denote the first and last equation of (2.6). We have the 
following equivalence of the two problems: 

Lemma 2.3 (Equivalence). When minr^x) > —1 the steady water-wave prob- 
lem (2.3) is equivalent to that £((r],^p),a,Q) = for 

M) G C 2 +£($,R) x {cJ+f CTen (A,M): ^| s =o = m ,^|.=i = m x } , 

and £((0, ip),a,Q) =0 /or ^ = ^(s) i/ and onfr/ if ip( s ) = V^ 5 )- 

Proof. Since 1 + 77(2:) 7^ 0, 77 is of class C 2+ " , and Holder spaces are Banach 
algebras, the identification ip(x,s) = ip(x,y/(l + rj(x))) = ijj(x,y) gives the same 
regularity of the functions ip and ip up to order 2 + /3. The equivalence then follows 
from direct calculation, and one can see from (2.6) that its x-independent solutions 
are also given by the formula (2.4) with y replaced by s. □ 

3. Functional- analytic setting. 

The map T . We want to linearize the problem £((r],ip),a,Q) = around a 
trivial solution tpg, whence we introduce a disturbance </> through ip = ipo + cf>, and the 
function space 

X := X 1 x X 2 := C^(§) x {4, G C*+f cvcn (^) : 4>\ s=1 = 0| s =o = o}. 
Furthermore, we define the target space 

Y := Y x x F 2 := (£J£(S) x <^ er> evei A 



and the sets 



and 



:= {(rj,4>) E X : minr? > -1} 



W := {(/x, a, A) € R 3 : \i ^ 0, a < 0, sin(A) ^ 0}, 
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which conveniently captures Assumptions 2.1 and 2.2. Then O C X is an open 
neighbourhood of the origin in X, and the embedding X K is compact. Elements 
of y will be written u> = (rj, 4>) and elements of IA will be written A = (/i, a, A). Define 
the operator T : O x U -> F by 

^(ifl.A) := (Ji( W) A),J- 2 (w,A)) 




sr ?x ((V'o) s (s;A) + ^ s ) \ ((^ ) s (s;A) + S ) 2 



1 + 7? 



(l + r/)2 



s=l 



and 



•F 2 (iu,A) 



s?7x((V'o) s (s; A) + 4> s 



S1]x 



0.i 



1+?/ / 

/ X 



I + 77 

a (tp {s; A) + 



1+7? 

(tA>Ms;A) + (^ 

(l + ?? )2 

We write C"^ to denote analyticity between Banach spaces. 

Lemma 3.1 (Equivalence). The problem J-((r), 4>), A) = 0, (77,0) £ 0, is equiva- 
lent to the water-wave problem (2.3). The map (r],(f>) <— > (77,1/') * s continuously differ- 
entiable, J"((0, 0), A) = «/ and oniy if(j> = 0, and T G C W (C x F). 

Proof. Relying on Lemma 2.3 it suffices to note that F((r], </>), A) = f ((77, t/>o(s; A)+ 
0), a, Q(A)). The analyticity follows from the fact that ip g C"(R x R), and that 
J 7 depends polynomially on <fi and ^o> whereas it is a rational function in 1 + n. □ 

Linearization. We obtain the linearized water-wave problem by taking the Frechet 
derivative of T at the point w = 0: 



D w J"i(0,A)u; = (V>o) s $ s - (tAo)s ? 7 + 77 



5=1 



(3.1) 



and 

D lu .F 2 (0, A) to = + SS -acf>- s(ip ) s ri xx - 2(ip ) ss r), (3.2) 

where D m denotes the Frechet derivative with respect to the first variable w. It turns 
out that there is a simple transformation that transforms this linearization to the one 
which is obtained by formally linearizing directly in (2.3). To that aim, let 

X 2 := [cf> e Cp^cvcn^): <£| s =o = o}, 

and 



X := {{ V J)eX 1 x X 2 } 
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Then X 2 X 2 and X X <-t Y, where the last embedding is compact. 

Proposition 3.2 (The T- isomorphism) . The bounded, linear operator T(A) : X 2 
X given by 

is an isomorphism. Define C(A) :— D w J-(0, A)T(A) : X 2 — > Y. Then 



C{K)<t>: 



OlsVs 



Wo) s 



(V>o) s 



, {d 2 x +d 2 s -a) 



(3.3) 



Proof. Note that we have the relation 

(pi + d 2 s - a)(s(tp ) s rj) = s(tp ) s r/ xx + 2(i/j ) ss r\. 
Thus, the right hand side of equation (3.2) can be written 

[dl + dl - a){4> - s(i> )s V) = {dl + d 2 s - a)4>, 
where cp = (p — s(ipo) s r\. Moreover, for <p <= X 2 we have that 

Y> w Fi{Q,N)w = [{ipo) s <p s + (i/'o) s (V'o) ss ?7 + 77]| s= i, 
<p\s=i = ~(ipo)s{l)v, 
<t>\s=o = 0. 

Due to Assumption 2.2 we can thus express r] in terms of <p and obtain that 

1 



T> w Ti(P,A)w 



(lpo)s <Ps ~ (ipo)ss + 



Wo), 



s=l 



REMARK 3.3. When the dependence on the parameters is unimportant, we shall 
for convenience refer to D w J^O, A), C(A) and T(A) simply as D w J-'iO), C and T. 

Via T, elements <p S X 2 can be "lifted" to elements (77, <p) £ X using the corre- 
spondence induced by the first component of T<p. The following result is immediate. 

Proposition 3.4 (Surface projection). The mapping r)r.) defined by 



X 2 3 tp h-» 7/4, 



W-o).(i) 

is linear and bounded, whence the linear map 

(?](.), id) : X 2 9 4> ^ (r?^,^) 6 X 

is bounded. 

We now have the functional- analytic setting outlined in Figure 3.1. Any property 
of the operator D^J^O) can be conveniently studied using the operator C. Some 
properties of the operator C are recorded in the following lemma. The proof is stan- 
dard and is therefore omitted (see e.g. [19, Section III.l] and [24, Lemma 3.5 and 
Lemma 3.6]). 

Lemma 3.5 (Fredholm properties). The operator C: X 2 — > Y is a Fredholm 
operator of index 0. Its kernel, ker£, is spanned by a finite number of functions of 
the form 



4>k(x,s) 



cos(kx)sm*(9 k s)/e k , 6 k ^0, 
cos(kx)s, 9 k = 0, 
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k € kZ, 



(3.4) 



1 



X C X C Ysfa*) 

T (J? ( .) ,id) 



£ = d u ,jf o r 

Fig. 3.1. The functional- analytic setting. 



where we have used the notation 



sin (9 k s) := 



sin(# fc s), fc 2 + a < 0, 
sinh(0fcs), fc 2 +a>0. 



Define Z := {{f](f>, (f>) : G kcr£} C X C Y. Then the range of C, ran£, is the 
orthogonal complement of Z inY with respect to the inner product 



(wi,W2)y ■= // <j>i4>2dxds+ / 7]ir] 2 dx, 
The projection Hz onto Z along ran£ is given by 



wi,w2 e Y. 



W,W k Y - 



Ffcllr 

where the sum ranges over all w k = {f](f> k , <frk) G Z, with (j) k of the form (3.4). 

The abbreviation sin* in Lemma 3.5 will be used analogously for other trigonomet- 
ric and hyperbolic functions. Noting that ranD^ .F(O) = ran£ and that kerD„, .F(O) = 
TkerZ! we obtain the following corollary of Lemma 3.5. 

Corollary 3.6. D w J-(0) : X — > Y is a Fredholm operator of index 0. 

Remark 3.7. If Assumption 2.2 is violated, that is if (ipo) s (l) = 0. then 
D tu J r i(0)u' = 7]. Consequently, the range of D w J- 1(0) equals X\, and this can be 
used to show that T) w J-(0) is not Fredholm. 

4. Bifurcation. This section contains the main proofs of the one- and two- 
dimensional bifurcation. 

Lemma 4.1 (Bifurcation condition). Let A = (//, a, A) G U, so that that Assump- 
tions 2.1-2.2 hold. For k G kL we have that cos(fcx) sm* (9 k s)/9 k G ker£ if and only 
if 



9 k cot*{9 k ) 



1 



yflBl sin (A) 



6» cot(A). 



(4.1) 



Remark 4.2. When 9 k =0 we interpret sm*(9 k s)/9 k as s and 9 k cot*(9 k ) as 1. 
We shall use this convention from now on. 

Proof. Assume first that 9 k ^ 0. By substituting </> = cos(kx) sm* (9 k s) /9 k into 
(3.3), we find that C(f> = if and only if 



9o/isin(A)cos*(0fc) = # /icos(A) 



1 



(i9q sin(A) 



.) 



(4.2) 



We claim that under Assumption 2.2 this equality implies that sin* (6^) 7^ 0. Indeed, 
if sin*(6>£.) = 0, then cos*(6*fc) 7^ 0. Examining the left hand side of equation (4.2) 
we find that 9Q(ism(\) = 0. But this contradicts Assumption 2.2. Dividing by 
9ofism(\) sm*(8k)/9k in (4.2) we obtain (4.1). When 9k — we substitute (j) = 
cos(kx)s into (3.3) and similarly obtain (4.1) with left hand side 1. □ 

Note that since we have restricted ourselves to even functions, it suffices to take 
k > 0. In order to find nontrivial waves we assume that k > 0. 

The bifurcation equation (4.1) is a transcendental equation involving the three 
parameters fx, a, A and the variable k. For fixed values of fi, a and A, the number of 
different solutions k G n1+ gives the dimension of ker 7^(0, A). A complete descrip- 
tion of the set of solutions and its dependence on the parameters is at the moment out 
of our reach. The next lemma shows, however, that one can choose the parameters 
such that there is precisely one solution k G kZ + . We also show that the parameters 
can be chosen so that there are precisely two solutions k%, G kZ + . It is an open 
question if there are values of the parameters for which there are three or more solu- 
tions, that is, for which the kernel of ^(0, A) is at least three-dimensional. Note that 
this question only depends on the left hand side of (4.1). For example, if for some 
a£l and a < there are exactly three solutions to the equation 9k cot*(9k) = a, we 
can always adjust A and [i to make the right hand side equal to a (this is clear from 
the proof of Lemma 4.3) and thus obtain an example of a three-dimensional kernel. 
The fact that the function k t— > 9k cot* (6k) is not monotone when a < is mainly 
what makes the analysis of the bifurcation equation difficult. On the other hand, it 
is exactly this property that allows for multi-dimensional kernels. 

Lemma 4.3 (Bifurcation kernels). Let ki, ki <G kZ+ and A G U. 

i) For a.e. a and any k\ there are [i and A such that (4.1) holds only for k = k\. 

ii) For any A with cot(A) < and fci,fc2 such that k\ > k\ + jft 2 , there are 
a G (— fcf, — k\ — ir 2 ) and fj, such that (4.1) holds for k = kx, k%, and for no 
other k > ki- 

Hi) For any ki, ki such that k\ > k 2 + 3ir 2 and k\ — k\ 7^ (2rj + l)7r 2 , for all n G Z, 
there are a < —k\, [i and A such that (4.1) holds at least for k = fci, /c2- 
The points a such that there exist fi and A for which (4.1) holds for more than one 
k G k7j + are isolated and thus a Lebesgue null set. 

Remark 4.4. For a as given by ii) and Hi) in Lemma 4-3, the background laminar 
flow (2.4) and its derivative (ipo) y (y; A) has at least one zero (stagnation occurs). For 
a as in Lemma i), the presence of a stagnation point for some y G (0, 1) is 
guaranteed if a < —w 2 . 

Note that in the cases ii) and Hi) we generally only assert that the kernel is at 
least two-dimensional. However, if we choose k\ = K in ii) we obtain that the kernel 
is precisely two-dimensional. In Example 5.3 it is shown that one can also obtain a 
two-dimensional kernel with ki > k. 

Proof. The right hand side of (4.1), call it r(/j.), satisfies 



if we choose A such that cot (A) < 0. We set t := \a\, and study the function 



ran(r(/i)) = (9 cot(A), 00) 2 R. 
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The functions h\ and h 2 both tend to 1 as t — k 2 vanishes, whence the function h is 
continuous away from the singularities t — k 2 = n 2 ir 2 , n g Z + . 
Now, let x := -J\t — k 2 \. Since 



d , , , cosh(x) sinh(x) — x 

— (x coth(x)) = , >0, 0<x£ 

dx K v " smh 2 (x) 

d , , .. cos(x) sin(x) — x 
— (xcot(x)) = \ > 
dx sin (x) 



— (x cot (a;)) = v , 2 v / < 0, < x <£ iri 



it follows that h is strictly increasing in A:, as well as strictly decreasing in t (which, of 
course, does not apply across the singularities). Moreover, if we let n be the maximal 
integer n > such that n 2 7r 2 < t, then 

J (v^cc-t y/t, oo) , t < 7T 2 , 
ran ( /i ) = <( „ 

t-fe 2 <t 1 M, t = 7T 2 , 

for n = 0, and 

2 



I (Vi cot v^, oo) , £<(n + l) 2 7r 
ran (ft) = < 
^7r=<t-fe2<t v IR, t=(n + l) 2 7r 2 



ran (/i) = R, n = 1, 2, . . . , n — 1, 

n 2 7r 2 <t-fc 2 <(ri+l) 2 7r 2 

ran f/i) = R 

t-fc 2 <77 2 

otherwise. 

To prove ii), consider the set / of |a| = t such that 

Sir 



<t-k{ < \—j and ran(/i(f;fci)) = (l,oo). (4.3) 

Due to the monotonicity with respect to k, the only possible fc > ki such that 
h(t;ki) = h(t;k) must belong to the interval k 2 > t, i.e. to the hyperbolic regime, 
and there exists at most one such k for a given t. In particular, if t fulfills (4.3) and 
k\ > k\ + jft 2 , then k\ — t > 0. Since, in that case, 

1 < h(t; k 2 ) < k 2 coth(fc 2 ), (4.4) 

while h(t;ki) spans (l,oo) as t spans /, we may choose |a| = t such that hi(t;k\) = 
h 2 (t; k 2 ). We then choose fj, such that h(t; kx) = r(/j). 

We now prove that the points a for which (4.1) holds simultaneously for at least 
two different k's are isolated. First say that 

f(t):=h 1 (t;k 1 )-h 1 (t;k 2 ), h < k 2 , (4.5) 
satisfies f(t ) = 0, for t > k 2 , j = 1,2. Differentiation with respect to t yields that 

cot U/i - kfj cot (yt - k%) 



/'(*) = 



2y/t - kx 2\Jt-k 2 



2 



icot 2 Ut-k 2 2 \ -\cot 2 I 
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which, using the fact that /(to) =0, at t = t can be rewritten as 

v _ (fej-M) ^ifa;fci)(feifa;fci)-i) 




A similar, although not identical, calculation for /ii(t;fci) and hi(t\k-2) yields the 
same conclusion when k\ < to < k 2 . When k\ < k\ = to a direct computation shows 
that /'(to) = —1/6, so that the conclusion remains true (here we have extended the 
definition of f(t) in the natural way). Since, in view of the monotonicity of h for 
k 2 — t > 0, there are finitely many fcj's for which one could have equality with k\ in 
(4.5). Therefore 



and a possibly smaller e exists such that h(t\ ki) ^ h(t; kj) for |t — to\ G (0,e), i ^ j. 
Hence, the points t for which /(t) = 0, for some ki ^ kj, are isolated. 

This argument also proves i): For a.e. t = |a|, we have that h(t;ki) ^ h(t;kj), 
kj fci.We then pick A such that 0$ cot(A) < h(t; fei) and the appropriate /j, such that 
h(t;ki) = r(fi). 

iii) Since we are interested in a < — fcf, we study t > k\. By the assumptions 
there exists n G Z + such that 



On the interval (fc| + n 2 7r 2 , fc 2 + (n + l) 2 7r 2 ) the function h(t; k%) spans E, while h(t; k\) 
is bounded. It follows that h(to; k\) = h(to; k%) for some to € (fe 2 + rt 2 7r 2 , fc| + (n + 
l) 2 7r 2 ). Hence, we may find (at least) one t = \a\ > k 2 such that h(t; ki) = h(t; A2) G 
K, and then choose A and [i appropriately, so that (4.1) holds for k = ki, fe- D 

Lyapunov-Schmidt reduction. Let A* denote a triple (/i*,a*,A*) such that 
(4.1) holds and suppose that the kernel is nontrivial, so that 



with 4>* = cos(kjs) sm*(9 kj s)/9k j and < k% < ■ ■ ■ < k n . Let w* = T(A*)0*. >From 
Lemma 3.5 it follows that Y = Z © ran£(A*). As in that lemma, we let Tlz be the 
corresponding projection onto Z parallel to ran£(A*). This decomposition induces 
similar decompositions X = Z © (ran£(A*) n X) and X = ker J- w (0, A*) © Xq, where 
X = ft(ran£(A*) n X) in which 



Applying the Lyapunov-Schmidt reduction [19, Thm 1.2.3] we obtain the following 
lemma. 

Lemma 4.5. There exist open neighborhoods M o/O in ker 7^(0, A*), M of in 
X andW of A* in R 3 , and a function tp G C°°(N X U',M), such that 

F(w,A)=Q for weM + M, A G W, 



inf {\h{t Q -k 3 ) - h{t -k x )\ : h{t o ;k J )^h(t o ;k 1 )}>S>0, 



(2n + 1)tt 2 < ki - k 2 < (2ra + 3)?r 2 . 



This implies that 



k\ + (n + 1) 2 tt 2 < k\ + n 2 TT 2 < k\ + (n + l) 2 ir 2 < k\ + (n + 2) 2 tt 2 . 



ker£(A*) = span{<^, </>*}, 
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if and only if w = w* + ip(w*, A) and w* = t\w\ + • • ■ + £ n U>* G M solves the finite- 
dimensional problem 

$(i,A)=0 /or ieV, AeW, (4.6) 

m which 

$(t,A) :=n z ^(«;*+V^K,A),A), 

and V := {( £ R" : fiiwj + ■ ■ ■ + G TV}. The function ip has the properties 

V>(0, A) = and D w ip(0, A) = 0. 

Bifurcation from a one-dimensional kernel. Theorem 4.6 (One-dimensional 
bifurcation) . Suppose that 

dimkcrDtu J"(0, A*) = 1, 

and let kerDy, .F(0, A*) = span{u>*}. There exists a C 00 -curve of small- amplitude 
nontrivial solutions {(w(t),Jl(t)): < \t\ < e} of 

T(w, fj,, a*, A*) = (4.7) 

in O x R, passing through (w(0), ~p(0)) — (0,fi*) with 

w(t) = tw* + 0{t 2 ) 

in O as t — > 0. In a neighborhood of (0, /i*) in O x R these are the only nontrivial 
solutions of (4.7). The corresponding surface profiles have one crest and one trough 
per period, and are strictly monotone between crest and trough. 

Proof. We apply the local bifurcation theorem with a one-dimensional kernel by 
Crandall and Rabinowitz [10]. We shall repeat the details here in order to clarify 
the proof in the case of a two-dimensional kernel below. By Lemma 4.5 the equation 
J-(w,A) = is locally equivalent to the finite-dimensional problem $(i, A) = where 
t G R. Note that <f>(t, A) = $i(i,A)w* where w* = (77^., 0*), 4>* = T^ 1 (A*)w*, and 
$1 is a real-valued function. Since $i(0,A) =0we find that <f>i(t,A) = W(t,A) in 
which 

*(t,A) = / (d t ®i)(zt,A)dz. 
Jo 

For f ^ the equations $1 = and '5 = are equivalent. We now prove that 
8^(0, A*) 7^ and apply the implicit function theorem to ty. We have that 

8^(0, A*) =d^ 1 (0,A*). 

Moreover, 

d^dMo, A*) = n z T(0, A*)w* 

(B 2 w ^(0 7 A*)w*,w*) Y 

Hence, d^(0,A*) 7^ if and only if (D^ J"(0, A*)w*, w*) y ^ 0. Recall that 
D„ J"(0,A) = £(A)r _1 (A). Thus, 

D 2 WM J-(0,A>* =D Al £(A*)^+£(A*)D AJ T- 1 (A*K. (4.8) 
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Since C = D w J-(0) o T the second term in the right hand side of (4.8) belongs to 
ran D.^ J"(0, A*), and we find that (D 2 wfl J"(0, A*)w*, w*) Y = (D p £(A*)(f)* ,w*) Y ■ A 
straightforward calculation shows that 

D M £(A>* 

■0 O sin(A*)cos*((9 fc ) - 



oCos(A*) 



(fi*) 2 Q sin(A* 



cos(fca;), 



so that 
(D M £(A* 



(V'o).(l) 



~*\ 

,W }y 

6>o sin(A*) cos*(6» fc ) + ( -0g cos(A*) + 



(fi*) 2 9 sin(A* 



sm 



sin 



(4.9) 



Using the fact that sin* (9 k )/0 k ^ (see the proof of Lemma 4.1) and rearranging, we 
obtain that the right hand side in (4.9) vanishes if and only if 



; cot*(0* 



1 



8q cot(A* 



which contradicts (4.1). □ 

Remark 4.7. It follows from the implicit function theorem and the bifurcation 
equation (4.1) that as a and A are varied in a neighborhood of (a* , A*) one obtains a 
whole two-dimensional family of bifurcating curves which depend smoothly on a and 
A. 

Bifurcation from a two-dimensional kernel. When the kernel is two-dimensional 
we can still apply Theorem 4.6 by restricting our attention to functions that are 2ir jki- 
periodic in the x-variable. Let X^ k2 ^ denote the subspace (subset) of such functions 
in any (open set in a) Banach space X . In particular, the restriction J 7 ^ 2 ) := J-\Q(k 2 ) 
is well-defined, and, in view of that k% > k\, we have that ker Du, J-"^ 2 -* (0, A*) = 
spanju;?;} := T(A*) spanjc^}. An application of Theorem 4.6 and Remark 4.7 to 
j?( k 2) yields that the set of nontrivial solutions of 



T(w, /it, a, A*) =0 
I 2 is locally given by the two-dimensional sheet 

< fc2 > = {(TO 2 (t,o) J 7Z 2 (t,a),a): ((,a)eV 2 ,(/0}, 



(4.10) 



where V2 is an open neighborhood of (0, a*) in M. 2 . When k-± is not an integer multiple 
of k\ the same argument gives a set t* 1 ) of nontrivial solutions of (4.10) in £)( fel ) x K 2 . 
It is however also possible to obtain solutions that are neither in X^ 1 ^ nor in X^ k2 \ 
as we will now describe. 

Theorem 4.8 (Two-dimensional bifurcation) . Suppose that 



Define 



dim ker D w J"(0,A*) = 2. 

a-.= e kl cot*{e kl ) = e k2 cot*(e k2 



as the left-hand side of (4.1), and assume that that a (jL {0, 1}, or that 9 k2 
which case a — 1). 



(4.11) 

= (in 
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i) If^/ki £ Z there exists a smooth sheet of small- amplitude nontrivial solu- 
tions 



mixed 



{(w(t 1 ,t 2 ),-p(t 1 ,t 2 ),a(ti,t 2 ))- 0<tj + t 2 2 <e 2 } 



of (4.10) m O x R 2 , passing through (w(0, 0), /Z(0, 0), a(0, 0)) = (0,/i*,a*) 

io(ii , * 2 ) = *iwif + h w 2 + 0{t\ + t 2 2 ) 

In a neighborhood of (0, //*, a*) iuOxI 2 i/«e union mlxed ij u (feO contains 
all nontrivial solutions of (4.10). 

ii) Let 5 > 0. // k 2 /k\ E Z £/iere exists a smooth sheet of small- amplitude 
nontrivial solutions 

mixed ._ ^{r,v),-p,(r,v),a{r,v))-. 0<r<e,6< \v\ < tt - 6} 

of (4.10) in O x M 2 , passing through (uT(0, v), Jl(0, v), a(0, v)) = (0,/i*,a*) 
with 

w(r 7 v) = r cos(v)w\ + r sm(v)w 2 + 0(r 2 ). 

In a neighborhood o/(0,/i*,a*) in O x R 2 the union f lxed U contains all 
nontrivial solutions of (4.10) such that 5 < \v\ < tt — 5, where r cos(v)wl + 
rsm(v)w 2 is the projection of w on ker J-'wiO, A*) parallel to Xq. 
Remark 4.9. When t\ ^= and t 2 ^ we find solutions which are neither 
2ir/ki nor 2tt /k 2 -periodic. A more precise result could be obtained by studying the 
Taylor polynomial of <f> at (0,A*) of sufficiently high order as in e.g. [18, 17, 22, 27]. 
The computations are in general quite involved. We settle for the mere existence of 
multimodal solutions here. 

Proof of Theorem J^.8. Define w* := (r^» ,</>*) £ X, j = 1,2, and recall that 
Z = span{w 1; w 2 }. If IIi<I> = and n 2 <I> = $ 2 «>2! where 111 and II 2 denote the 

projections onto span{w^} and spanjw^}, respectively, equation (4.6) takes the form 

$i(ii,i 2 ,A) = 0, 

H 1» 2, , 

* 2 (ti,t2,A) = 0. 

This is a system of two equations with five unknowns, and we note that it has the 
trivial solution (0, 0, A) for all A g W . 

Assume first that k 2 /k\ ^ Z and suppose that t\ = and w* = t 2 w 2 , t 2 € R. 
An application of the Lyapunov-Schmidt reduction in the subspace of 2-7r/fc 2 -periodic 
functions yields that ip(w*,A) is then 27r/fc 2 -periodic. Hence, 

$i(0,i 2 ,A) = 0, for all t 2 , A. (4.13) 

Similarly, 

$2^1,0, A) = 0, for all ix, A. (4.14) 

Let * := (*i,* 2 ), 

*i(ti,ta,A):= / (d tl $i)(zt u t 2 ,A)dz, ^ 2 (h,t 2 ,A) := f {d t2 $ 2 )(h, zt 2 , A) dz. 



(i 
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From (4.13) and (4.14) it follows that (4.12) is equivalent to that 

ti*i(ii,i 2 , A) = 0, 

\ (4-15) 
t 2 *a(*i,*a,A) = 0. 

There are four possibilities. The trivial case t\, ti =0 corresponds to trivial solutions. 
When t\ — but t% 7^ the system reduces to ^(O, ti, A) = 0, the solutions of which 
are locally given by the set ( fc2 ); when t\ ^ but t2 = we similarly obtain the set 
( fcl ). It remains to investigate the solutions of \&i(ti, t%, A) = ^2(^1,^2, A) = in a 
neighborhood of (0,0, A*). Clearly, when those exist, they will intersect and connect 
the sets ( fcl ) and Since 



we have that 



Similarly, 



3 tl $(0,0,A*) = Il z D w F(0,A*)w* 1 =0, 
* 1 (0,0,A*) = 9 4l $ 1 (0,0 J A*) = 0, 
* 2 (0,0,A*) =3 t2 $ 2 (0,0,A*) =0. 



We now apply the implicit function theorem to at the point (0,0, A*). It suffices 
to prove that the matrix 

9^i(0,0, A*) 9^2(0,0, A*) 

a Q ^!(o,o,A*) a Q * 2 (o,o,A*) 

is invertible. To this effect, we note that 

9^x(0,0,A*) = 5 tl ^$i(0,0,A*), 



and 



We have that 



d tl d^(0, 0, A*) = U z F(0, A*)w* v 



9 (DL JY0,A*)iy?,w?W 

z D «,/(0> A )™i = ~ jnrjrp ™i> 

J i lly 



since J-"(0, A*)«;* is proportional to cos(fcicc) and therefore orthogonal to w^- 
Consequently, 

(D, 2 „„ T(Q,A*)w?,w?) Y 
a M *i(0,0, A*) = 3^*1(0,0, A*) = ^ ^ V / 2 - 17 ■ 

ll^i lly 

Using the observation (4.8) we find that 

W0,0,A) = 

Similar arguments show that 

^ 2 (o,o,A)- , 
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and 



d Q ^(0,0,A*) = 



3 = 1,2. 



Thus 

dot 



M *i(O,O, A*) 5^*2(0,0, A*) 
a *i(O,O,A*) a Q * 2 (0,0,A*) 



= Cdet 



(D Al £(A*)^,^) y (D M £(A*)^,tS5>y 

(D a £(A*)#, ti^y (D Q £(A*)0*., W 2 *)y 



where C = ||to* ||y 2 ||w2lly 2 7^ 0- From (4.9) we see that 
(D„£(A*)0*,^) y 

9 sin(A) cos* (6 k ) + ( -0 2 cos(A) + 



(^b).(l) 

7T 



(/z*) 2 o sin(A) 



sin (6^, 



?o sin(A)a — 6q cos(A) + 



1 



( M *) 2 #osin(A) 



A 



A straightforward calculation shows that 

fj, sin(A) 



Da£(A*)0 = 

From this it follows that 



2fln 



4> s — I (j, cos(A) 



2/isin(A)6»g 



<D a £(A*)^?) y =iJ^^) + /(%), 



where _B 



W>o)»(i) 



^ sin A 

29n 



a — ( \i cos A — 



2fi sin A0q 



and 



/(%) 



7T /, 2 i -cos*(e fc )sin*(6» fc ) 

f sgn(fc| + a) — 2 p 2 L , 6> fcj ^ 0, 



sin*(0 fcj ) 



We thus have that 
det 
= det 
= A det 



A 



(D P £(A*)</,*,^*) F (D M £(A*)^,«)5) y 
(D a C(A*mwl) Y (D a £(A*)^,* 2 *)y 

A(sin*(0 fel )) 2 /^ i A(sin*(0 fe2 )) 2 /^ 2 
BCsin*^)) 2 /^ B(sin*(^))7C + /(fc 2 ) 

(sin*^)) 2 /^ (sin*(0 fc2 )) 2 /0 2 2 
f(ki) f(k 2 ) 

(sm*(9 kl )f ffu , (sin^)) 2 



? 2 
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Note that if a = 0, then fc 2 + a < and cos(0fc 3 .) = 0, j = 1,2. Hence, f(kj) = 
—ir/(28l ) in that case, and 

- ^p^-m) = - * - S i^(^)) = o. 

k\ k-2 k\ ki 

We therefore assume that a ^ from now on. We consider the three different possible 
cases. 

Case I. Assume that kf + a < and k\ + a < 0. In this case we find that 

^2 /( fc 2) ^ 

7T 1 



2^ 2 



a" 1 [sin 2 (0 fe2 ) ((a - 1) + sin 2 (# fcl )) - sin 2 (0 fel ) ((a - 1) + sin 2 (0 fca ))] 



f ^V a " 1(a " 1} (sin 2 (e fe2 ) -sin 2 (0 fel )) 

k\ fc 2 



This is non-zero if a ^ 1, in view of the equality 0^ cot^^) = #fc 2 cot(0fc 2 ). 

Case II. Assume that k\ + a < and k\ + a = 0. In this case a = 1 and we find that 

(sin*(0 fel )) 2 . (sin*(0 fc2 )) 2 ^sin 2 (0 fel ) tt / cos(0 fcl ) sin(0 fcl ) 



_ tt sin 2 (6> fcl ) 

~ 6 °l ' 

fci 

which is non-zero by the argument used in Lemma 4.1. 

Case III. Assume that k\ + a < while k\ + a > 0. We then find that 

(sin*(fl fcl )) 2 (sin*(g fc2 )) 2 

% °k 2 
TT 1 



a sin 



fel ) ((a - 1) - sinh 2 (0 fc2 )) + sinh 2 (0 fc2 ) ((a - 1) + sin 2 (0 fel ))] 



2 2 0? 

= |fl?V - 1) (sin 2 (^ fel ) + sinh 2 (^ 2 )) ^ 0, 

fei fc 2 

when a 7^ 1 . This concludes the proof of the first part of the theorem. 

Assume next that kilk\ € Z. In this case (4.13) remains true, whereas (4.14) 
may be false. We introduce $i as before, but write it in polar coordinates: 

$?i(r,v,A) = J 9i 1 $i(zrcos(u),rsin(u), A) dz. 
Jo 

We do not know that $2 vanishes when t% = 0, but we can still use the identity 
$2(0,0, A) = 0. Redefine * 2 as 

* 2 (r,u,A) 

^dt 1 §2(zr cos(u), zr sin(u), A) cos(w) + dt 2 $2(zr cos(u), zr sin(u), A) sin(u) j dz. 
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We then have that r^ 2 { r iV,^-) — &2{ r cos(u), r sin(w), A), and equation (4.12) is 
equivalent to (* 1 (r, u, A), \E , 2 (r, u, A)) = (0, 0) whenever t 1 ^0. Moreover, *i(0, v, A) = 
^2(0, v, A) = for all A. The complete set of solutions is locally given by the three 
cases r = 0, cos(w) = with ^2(r,v, A) = 0, and Vt , i(f, v,A) = ^(r, v, A) = 0; 
the first set is trivial, the second given by ^' 2 ', and the third yet to be found. By 
considering <9 r $(0,v, A*), as before we have that 



Filly 



while 



OA-,,,, » *\ (D M £(A*)#,rii 2 )y (D M £(A*)^,ri);)y . 

fy*a(0, w, A ) = — -p cos(u) + sin(u) 

II <"2 II y II "'ally 

where we have used the fact that D M £(A*)</>* is orthogonal to u>|. Similarly (cf. 3.3), 

Q T ,„ a *"\ (D Q £(A*)</>*,™*)y j _ . ... (D Q £(A*)0*,iZ;*)y . , , 
Q *i(O,u,A*) = ^ , .;.; 2 17 and a Q * 2 (0,u,A*) = ^ , ,;.; 2 2 ' 2/ am(v). 



j i\\y n^2iiy 



Thus, 



d^ 1 (o,v,A*) 0^2(0, u, A*) 
a Q *!(o,u,A*) a Q *2(o,w,A*) 



det 

= C sin(w) det 



(D M £(A*)#,u;{>y (D M £(A*)^ 2 ,«)5>y 

(D Q £(A*)<^,^)y (D Q £(A*)0 2 ,TD5)y 



where C = ||w*|jy 2 ||?5) 2 |jy 2 7^ 0. We can therefore apply the implicit function theorem 
under the assumptions of Theorem 4.8 as long as sin(v) 7^ 0. This gives us local, 
smoothly parametrized, solution sets at every point (r, v) = (0, vq) with < vq < it. 
On each compact interval in [S, Tt — S] C (0, tt) we may find a positive number e such 
the parametrization is uniformly valid for < r < £, 5 < v < tt — 5. □ 

Remark 4.10. By varying X in a neighborhood of X* we obtain a smooth family 
of bifurcating two-dimensional sheets of non-trivial solutions. This is quite natural if 
one considers the bifurcation equation (4.1): if X is varied near X* one can adjust fx to 
make sure that the kernel remains two-dimensional. One can check that the theorem 
also holds if (A, a) are used as bifurcation parameters instead of a). However, the 
vectors (<9 M *i(0, 0, A*), d^ 2 (0, 0, A*)) and (d x ^i(0, 0, A*), A *a(0, 0, A*)) are paral- 
lel, and it is therefore not possible to replace {is, a) with (n,X) (the same problem 
appears when k^/ki G 7L). Again this can be explained by looking at (4.1). If the 
result were true with (//, A) as bifurcation parameters, we would obtain a family of bi- 
furcating two-dimensional sheets indexed by a. However, Lemma 4-3 shows that if a is 
varied near a* , it is impossible to adjust \i and X to keep the kernel two-dimensional. 

5. Applications and examples. The waves found in this paper are all small- 
amplitude rotational gravity waves in water of finite depth. For simplicity we shall 
take k = 1 here so that the basic period is 2ir. 

Corollary 5.1 (Stokes waves with arbitrarily many critical layers). Let j3 <G 
(0, 1), k £ Z + . For almost every a < —tt 2 there exists a C°° -curve 



t 1 V 



ri,i>,Q,m ,mi] (t) e C 2+/3 (§) x C 2+/3 (fi) x 
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Fig. 5.1. Rotational gravity water waves bifurcating from one- and two- 
dimensional kernels. Left: the qualitative features of surface profiles for the small- 
amplitude waves found in Corollary 5.2 along the bifurcation parameter v. The marked lines 
show: at v = tv/2 a typical Stokes wave with a single crest within each period; at v ~ 2n/3 a 
wave with several crests within each minimal period. Right: one typical rotational behavior 
of the waves found in Corollary 5.1 bifurcating from a single eigenvalue. The dotted lines 
indicate streamlines. Notice that the number of vortices can be chosen arbitrarily large. 



corresponding to even, 2ir/k -periodic solutions (ry, ip, Q, tuq, mi) of the steady water- 
wave problem (2.3). These have a surface profile which rises and falls exactly once in 
every minimal period, and at least one critical layer. The number of critical layers 
can be chosen arbitrarily large by choosing —a large enough. 

Proof. Lemma 4.3 i) yields the existence of (i* , A* and, via Lemma 3.5, the one- 
dimensional kernel w* — T(A*) cos(fca;) sin* (8ks)/9k needed for Theorem 4.6 (here 
a* := a). By identifying m\ := fi(t) cos(A) and m.2 := ^(t)cos(A — ^/ja*]), we 
obtain from Lemma 3.1 the solutions of the original problem, as well as the (In- 
dependence with respect to t. Notice that also Q = Q(/i(t), a*, A*) depends on the 
bifurcation parameter t through the formula (2.5). If we identify w(t) = T{t)4>(t) we 
have D t n(0) = r)(p t ^ )y. 



D 



( -m 

Wo)s{t) 



t=0 



f Vtttt) 1 



= D t (HO) 
t=o W<o) s (0) : 



since t = corresponds to the trivial solution (f>(0) = 0. In view of that 

sin*(6> fc )cos(fcx) 

we have C(A*, k) ^ and the local expression 

r)(x;t) = tC'(A*,k)cos(kx) + O (t 2 ) , as t ->• 0, 

in the C 2+ "(E>, K)-metric. This guarantees that the surface profile has one crest and 
one trough per period for t small enough. The number of critical layers for small 
enough solutions along the bifurcation curve can be seen directly from the number of 
zeros of (ipo)y = — M^o sin(6*o(y — 1) + A) (for details, see [11]). □ 

Corollary 5.2 (Doubly-periodic waves with at least one critical layer). Let 
(3 £E (0, 1) and k > 5 be an integer. There exists a family 

(r,v) i — ^ [77,V,Q,m ,mi,a] (r,v) € C 2+/3 (§) x C 2+/3 (H) x R 3 x (-fc 2 ,-7r 2 - 1), 

of nontrivial, even, 2it -periodic solutions of the steady water-wave problem (2.3) with 
at least one critical layer. For each v £ (0,7r) their surface elevation satisfies 



rj(x; r) = rC\ cos(u) cos(x) + rCi sm.(v) cos(kx) + O (r 2 ) , Ci,C2^0, (5.1) 
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in the C 2+p (S,R) -metric as r -> 0. 

Proof. The assumption on k guarantees that we may pick any A* with cot(A*) < 
and use Lemma 4.3 ii) with k\ = 1, k 2 = k, yielding fi* and a*. The proof of 
Lemma 4.3 ii) is devised so as to have a £ {0, 1} (this is (4.4)). That, in turn, allows 
us to apply Theorem 4.8 ii). The constants mo and m% are identified as in the proof 
of Corollary 5.1 where now also a(r,v) depends on the bifurcation parameters, and 
(2.5) means that Q = Q(jl(r, v), a(r, v), A). Since a* e (— fc 2 ,— ir 2 — 1) according to 
Lemma 4.3, and 5(0, v) = a* for all v £ (0, 71"), we still have a(r, v) G (— k 2 , —ir 2 — 1) 
if, for each v € (0, it), we choose \r\ small enough. The formula (5.1) is obtained from 
d r w(0, v) — cos(i;)T(A*)(/>* +sin(w)T(A*)</>2 just as in the proof of Corollary 5.1, and 
the existence of a critical layer follows from (2.4) and the fact that a < —ir 2 . □ 

Example 5.3 (Doubly-periodic waves with several critical layers). There exists 
a family of nontrivial, even, 2tt -periodic solutions of the steady water-wave problem 
(2.3), with approximate surface elevation 

r)(x) = cos(4ir) + C 2 t 2 cos(7x), C x , C 2 ^ 0. 

and three critical layers. 

Sketch of proof. One checks that k\ = 4, k 2 = 7, and a* ~ —60 fulfills (4.11) for 
a w 18. Those k^s are thus an example when Lemma 4.3 hi) holds exactly for k\ and 
k 2 . The proof then follows the lines of that of Corollary 5.2, except that we now use 
Theorem 4.8 i). Since yj\a* \ > 2tt it follows from (2.4) that, for t 2 +t 2 small enough, 
the resulting waves have at least two critical layers. By choosing first A and then /i 
appropriately in (4.1) one obtains three critical layers. □ 
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